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ABSTRACT We define the osp(l,2) Gaudin algebra and consider integrable 
CT) | models described by it. The models include the osp(l,2) Gaudin magnet 

and the Dicke model related to it. Detailed discussion of the simplest cases 
of these models is presented. The effect of the presence of fermions on the 
separation of variables is indicated. 

<D 

> : 

. I. INTRODUCTION 

X 

In 1976 |l| (see also 0) M. Gaudin introduced a quantum mechanical model of N particles 
with su(2) or su(l, 1) spin, ever since known as a Gaudin model. The model is described 
by a set of iV hamiltonians 

Hi = ^=-^-, i=l,...,N, (1) 

where e\ > e<i > . . . > are free parameters and S l is the su{2) or su(l, 1) spin of the 
z-th particle, i.e. 

[Si, Si}- = ±6ijSi, [S^ = 25 tJ Sl. (2) 

Here [ , ]_ denotes the commutator. The total spin S 2 = I]j(S 4 ) 2 and S z = J2iSl are 
remaining constants of motion. Gaudin showed that the hamiltonians (|I|) commute with 



each other and constructed their common eigenvectors by the coordinate Bethe Ansatz. 
In || the Gaudin model was generalised to any semisimple Lie algebra. 

Because of its simplicity, the Gaudin model has been used as a testing ground for 
the ideas such as the functional Bethe Ansatz and the general procedure of separation 
of variables ||]. Its connection to the Neumann model ||, revealed in || allows one to 
generalise separation of variables for some types of coordinates on Riemannian manifolds. 
It has also applications in quantum optics [^]. 

In [|J the Gaudin model was solved in the framework of the algebraic Bethe Ansatz. 
It was also shown there [|J that the model is governed by a Yang-Baxter algebra, called a 
Gaudin algebra, with commutation relations linear in the generators and determined by 
a classical r- matrix. It is to be stressed these features are present in the model despite its 
quantum mechanical nature. In fact, the Gaudin model is one of a large class of models, 
with such an algebraic nature, so that study of Gaudin algebras becomes an important 
issue. 

To construct a Gaudin algebra we need two ingredients: the classical r-matrix and 
quantum L-operator. By the former we mean an M 2 x M 2 matrix, of which the compo- 
nents are functions satisfying the classical Yang-Baxter equation 

[n 2 (A - /i), n 3 (A - u)]- + [n 2 (A - //), r 23 (/i - z/)]_ + [r 13 (X - u), r 23 (> -!/)]_ = 0. (3) 

The subscripts in (^) indicate the vector spaces in R M ® R M ® R A/ on which r acts 
non-trivially, and A, /x, v are spectral parameters. The quantum L-operator is an M x M 
matrix, of which the components are operator valued functions, such that 

[L x (A), L 2 (/i)]_ = -[r(A - n), la (A) + L 2 (/i)]_, (4) 

where Li(X) = L(A)£g>l and L 2 (/x) = l®L(fi). The infinite dimensional algebra generated 
by the components of L(X) is called a Gaudin algebra. If a Gaudin algebra describes the 
Gaudin system corresponding to Lie algebra Q, then we term this kind of a Gaudin 
algebra a Q Gaudin algebra. To any simple Lie algebra one can associate different types 
of Gaudin algebras, corresponding to different types of r- matrices [Q. For instance we can 
consider rational, trigonometric or elliptic Gaudin algebras. In what follows we restrict 
ourselves to the rational case. Using language just introduced we can say that the model 
([!]) is described by the rational su(2) or su(l, 1) Gaudin algebra. 
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It is easy to extend the notion of a Gaudin algebra to superalgebras. In this case, 
both r(A) and L(X) have a supermatrix structure however. Hence the tensor product in 
(|3|) and (|4|) must be understood in a Z 2 -graded sense. Precisely, Eq.(Q) can be written as 

L(X) ® L{ji) - (-l) 9i W 9l ML(/i) <g> L(A) 

= L(A)r (1) (A - /i) (8) r (2) (A - //) + (-l) 9L ^*«r (1) (A - //) ® L(//)r (2) (A - /i) 
- (-l) aL(A)a ^r (1) (A - ^)L(A) ® r (2) (A - /x) - r (1) (A - //) ® r (2) (A - /i)L(/i), 

where we have used the notation r(A — n) = T(i)(A — ji) ® T(2)(A — /u) and e.g. dL(fi) 
denotes the parity of an appropriate element of L(n). 

In the present paper we describe models which can be solved using osp(l, 2) Gaudin 
algebra. In particular we generalise the Gaudin model by considering a system of N 
hamiltonians 

H< = £ »-y + iq*-vM j = li ... iW> (5) 

where the components of S*, satisfy the algebra @ and 

^]+ = ^^4. tn> ^]+ = ^^5j, (6) 

where [ , ]+ denotes the anticommutator. 

The paper is organised as follows. In Section II we construct the osp(l, 2) Gaudin 
algebra and we define a generating function for integrals of motion of systems described 
by this algebra. In Section III we construct the spectrum of the generating function 
and we derive the Bethe equations. In Section IV we use the osp(l, 2) Gaudin algebra to 
solve the osp(l, 2) Gaudin model. We show that the hamiltonians (^]) commute with each 
other and construct explicitly their complete spectrum, treating the two particle (N=2) 
special case as an example. In section five we discuss the osp(l,2) Dicke model. This 
is a natural osp(l, 2) generalisation of the ordinary Dicke model governed by an su(l, 1) 
Gaudin algebra, a generalisation which is realised by coupling a system of boson and 
fermion oscillators to the Gaudin model itself. We prove that the spectrum obtained for 
it by the algebraic Bethe Ansatz of section three is complete. Then we focus our attention 
on the special (N=l) case in which there is only one particle in the underlying Gaudin 
model. In an oscillator representation of osp(l,2), this case describes a system with 
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two oscillators with frequency ratio 2:1 together with one ordinary and one Majorana 
fermionic oscillators. We discuss the non-separability of the Schrodinger equation for this 
case. 



II. STRUCTURE OF THE osp{l,2) GAUDIN ALGEBRA 

Here we define the osp(l, 2) Gaudin algebra. To do it we first define a classical r-matrix 
and a quantum L-operator. The r-matrix is constructed out of the quadratic Casimir of 
osp(l,2) in a standard way ||, 



r(A) = j(S z 



S Z + -{S + ®S^ + S^ 



s+) + v+ <g> f_ - y_ ® v+ 



A G R, 



(7) 



written in the fundamental representation of osp(l, 2) (unitary in the compact case, non- 
unitary in the non-compact case), 
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Here and in what follows we use the convention, in which a matrix element is even 
(odd) if i + j is even (odd). This convention proves more suitable for our purposes, than 
the standard one, in which a supermatrix is written in a canonical block form (cf. ||10|| ). 
Inserting (H) into (I7|) we obtain an explicit expression for r(A) 



''< A > = 2A 
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Next we define the quantum L-operator 

/ A(X) V4X) B(X) \ 
L(X) = -V+(X) V-(X) . (10) 
V C(X) V+{\) -A(X)J 

The components of L(X) generate the osp(l, 2) Gaudin algebra. We assume that they are 
non-singular except for a finite number of distinct values of A, namely e\, . . . , ejy, such 
that ei > . . . > cat. 

Using definitions @, fllQl) and relation <Q) we obtain the complete set of relations 
defining the (rational) osp(l,2) Gaudin algebra 

[A(X), A(//)]_ = [B(A),B(ai)]_ = [C(A),C(//)]_ = 0, 

[B(X), V-((i)]-. = [C(X), V+(fi)]- = 0, 
(7(A) - C{p) 



[A(A), £(//)]. 
[A(A), £(//)]. 



A — 

fi(A) - B(ai) 



[A(A),V±(a*)]_ = =F 



A — /x 

iy±(A)-y±(/i) 



2 A - n 
■(A)-Afc) 

X — n 
1 A(X) - A(jjl) 
' 2 A - /i ' 
V_(A)-V-(/x) 

A — /j 
y + (A)-y + (/i) 
A — /! 
1 £(A) -£(//) 
2 A-// 
1 C(A) - C(/i) 



[B(A),C W ]_ = 2^iM, (11! 

[K_(A),y + (Ai)]. 
[s(A),y + 0i)]_ 

[C(A),V_(a*)]_ 
[VL(A),y_(Ai)]. 

[F + (A),y + (/i)]. 



2 A-// 

In the osp(l,2) Lie superalgebra there is a clear relation between V± and S±, i.e. 
4V± = T'S'i- Therefore the universal enveloping algebra C/(osp(l,2)) is generated by V± 
and 5* 2 only and the highest weight representations are constructed by the actions of a 
single raising operator V + say. In the Gaudin algebra the similar relation can be read off 
the last two of Eqs. (|TT|). When /i approaches A we obtain 

B'{X) = 4K_ 2 (A), C\X) = -4\/ + 2 (A), (12) 
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where prime denotes differentiation with respect to A. We learn from Eqs. ([12] ) that 
construction of highest weight modules of osp(l,2) Gaudin algebra involves both V+(A) 
and C(A) operators. This also makes the Bethe Ansatz construction of the next section 
more complicated, as compared to the su(l, 1) case, in which there is only one raising 
operator. 

From the point of view of integrable models described by the osp(l, 2) Gaudin algebra, 
it is important to consider a generating function 

t(X) = istrL 2 (A) 

= A 2 (A) + ^5(A),C(A)] + + [F + (A),\/_(A)]_. (13) 

This definition implies 

[t(A),*0u)]_ = O, (14) 

so that the function £(A) generates a family of commuting hamiltonians. Solution of an 
integrable model is equivalent to the construction for it of the spectrum of £(A). This 
can be carried out in the framework of the algebraic Bethe Ansatz, as will be discussed 
in the next section. 



III. THE SPECTRUM OF t(A) 

Now we construct the spectrum of the generating function t(X). We look for the highest 



weight representations of the Gaudin algebra (|il|) with a vacuum |0) defined by 



V1(A)|0) = J B(A)|0) = 0, A(A)|0) = a(A)|0). (15) 
Using definition (p~3|) of t(X) it is easy to find that 

t(A)|0) = (a 2 (A) + i«'(A))|0) 

= T (A)|0>, (16) 

where the prime denotes differentiation with respect to A. Other eigenstates of t(X) 
are then generated from the vacuum |0) by the repeated actions of V+(A) and C(A). 
Accordingly we look for an eigenstate . . . , fi n ) of £(A) of the form 

\fl 1: . . . , fl n ) = 1/ + (/ii,...,/i n )|0). (17) 
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Here V + (fi\, . . . , fi n ) is constructed from V+(//j) and C(//i), z = l,...,n. Demand- 
ing that \fjLi, . . . , fi n ) be an eigenstate of t(A) we determine the form of V + (/ii, . . . ,/i n ) 
and the numbers /ii, . . . ,/x n . The corresponding eigenvalue of t(A) will be denoted by 
r(A; /ii, . . . , /i n ). The explicit construction of V + (fii, . . . , /z n ) is rather involved. We list 
operators V + (//i, . . . , /i n ) for the first few excitations at the end of the section. Neverthe- 
less it is possible to find the eigenvalues r(A; /ii, . . . , /i n ) as well as numbers fi\, . . . , \i n 
without performing a fully explicit construction of V + (//i, . . . , fi n ). 
The operator V + (fii, . . . , /i n ) may be decomposed as follows 

V+(/ii, . . . , /i n ) = VVfjUl) • • • V^fjUn) + . . . , //„), (18) 

where W(/xi, . . . ,/i n ) involves at least one operator C(//j), i.e. 

1^0*1,...,^)= 51 7 CT (/^l,---,^)^+(^<Tl)---^+(^„- 2 )C , (^n-l) + --- ; 

a 1 <...<<7 n _ 2 

where 7o-(/xi, . . . , /i n ) are number coefficients. It is crucial to observe that W(fii, . . . , fi n ) 
cannot produce terms proportional to V+(fii) • • • V+(fi n ), when commuted with t(X). Like- 
wise W(jii, . . . , fi n ) itself cannot generate an eigenvector of t(X). Since ( |17D defines an 
eigenvector of t(X), the eigenvalue r(A; /ii, . . . , //„) has to be equal to the coefficient of 
V+(pi) ■ ■ ■ V+(/x n ), when acted upon by t(X). We compute 

n 1 1 1 

fziX-fii 2f^ j (X-n i )(X-fi j ) 

+ (terms involving at least one C(/Xj)) • (19) 

From this we immediately see that 

n 1 1 1 
r(A; . . . , /v) = £ t «(A) + - £ 7- + r (A). (20) 

The numbers //j, i = 1, . . . , n are not yet specified. They can be determined as follows 
[TT[| . We observe that the operator t(X) is non-singular for any X ^ €{, i = 1, . . . , N. In 



particular, fa, i = 1, . . . ,n are regular points of t(X). This implies that r(A; fix, . . . , fi n ) 
should be non-singular at A = fii, i — 1, . . . , n. The conditions for non- singularity in (p0|) 
at A = /ii, z = 1, . . . , n read 

1 n 1 

«W = "nE — Z — ' i = l,...,n. (21) 

2 j=i /^i A*j 
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These are the Bethe equations for the osp(l, 2) Gaudin algebra. Taking (|21|) into account 



we can write fl20|) as 



r(A; //!,..., /i n ) = r (A) + ± a(A > a(/Xi) . (22) 



We can also write (20) in the form 



r(A; /ii, . . . , // n ) = (x(A; ^i, . . . , /i„) + a(A)) 2 + \j^(x(^ Hi, ■ ■ ■ , Hn) + <*(A)), (23) 



where 



x(A;/ii, . . . ,/i„) = ttMtv, ?n(A) = II(A - fa). 

Z Qn{^) i= i 



Eq. (p3|) , rearranged as a differential equation for g n (A), reads 

g£(A) + 4a(A)g^(A) + 4(r (A) - r(A; . . . , /x„))g„(A) = 0, (24) 

so that (|20| ) and ( pi]) are equivalent to the differential equation (|24]) . This equation plays 
an important role in the analysis of the completeness of the spectrum of t(X) constructed 
by the Bethe Ansatz. We illustrate this in the next section. 

We would like to end this section with a few comments on the structure of operators 
V+(/j,i, . . . , /!„). Since each Hi corresponds to a fermionic degree of freedom, the function 
V + (hi, . . . ,/!„) is totally antisymmetric. The first three operators V + (hi, . . . ,/i n ) come 



out as 



4 Hi — H2 



v + Mv + Mv + M - 1 E e ijk v + ( tM ) C(fij) + C{flk) . 

8 i,j,k=i V-i ~ V* 

Here V+(hi) ■ ■ ■ V+(nn) denotes the antisymmetric product 

cr£S„ 

It can be checked directly that these operators generate eigenvectors of t(X), provided 
that conditions fl2~T|) are satisfied. 

Finally we notice that if there is no pseudovacuum fljjD, but there is a pseudovacuum 
V+(X)\0) = C(A)|0) = 0, A(X)\0) = 5(A) |0) instead, then the whole of the construction 
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can be repeated with V" + (A) replaced by — VL(A), C(A) replaced by —B(X) and a(A) 
replaced by —a (A). 



IV. THE osp(l,2) GAUDIN MODEL 

Here we discuss the model defined by the hamiltonians (^j). We show that hamiltonians 
([5]) are in involution and construct their common spectrum. We also prove, this spectrum 
is complete. Finally we derive the spectrum explicitly for the simplest non-trivial case of 
the two-particle system. We begin with a proof that the model is integrable. 
We consider the L-operator 



v 1 



si VI -SI 



Using the commutation rules (|2|) and (^), we can directly check that L(X) generates 
the osp(l,2) Gaudin algebra. Writing generators of the osp(l, 2) Gaudin algebra in the 
explicit form Q2"5D, we easily find that the generating function (|IBD takes the following 
form 

,m-V ( S ') 2 . ^ S< ■ Si + VjVl - VlVj 
t{X) -k^W + h (A-e,)(A- fJ ) (26) 
Hamiltonians (|5]) are identified with 

Hi = res A=e i(A) (27) 

and thanks to Eq.(|14|) they commute with each other. They also satisfy Y^iHi = 0, so 
that ( p7[) gives N — 1 independent hamiltonians. It is obvious that S 2 commutes with 
all Hi. Now it remains to show that S z commutes with t(X) and hence with each of H^s. 
The easiest way to do it is first to observe that the algebra ( |TT| ) will not change if we 
transform 

A(X)^A(X)+g, geR. 
This transformation induces the transformation 

t(X) -> t(A) +g 2 + 2A(X) = t(A). 

Obviously t(X) generates a family of commuting hamiltonians and 

Hi = res x= J(X) = H t + 2gS\. 
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Now J2i Hi = 2gS z and from [Hi, Hj]~ = we get [S z , Hi]- = as claimed. 

From now on we focus on the non-compact case, in which osp(l, 2) contains su(l, 1) 
as a subalgebra. We show, how to construct the complete spectrum of t(X) (^6|). We 
define the numbers s» > 0, i — 1, . . . , N by 

Si\0) = Si \0), i = l,...,N. 

With these definitions 

N 8- 

«(a) = Et^V 

i=l A ~ € i 

and the Bethe equations ( [Zip take the form 



N s- 1 " 1 

i=l A — C; 2 j=1 fii — flj 



n. 



They can be identified with those corresponding to su(l, 1) Gaudin model by the rescaling 
Si — ► Therefore, as in the sit(l, 1) case, the differential equation fl24|) becomes the 
Lame equation, intensively studied in the context of orthogonal polynomials [l^j. The 



result, which we need here is the Heine-Stieltjes theorem [0 asserting that there are 

r» 



'JV-1 



n + N -2' 

n 



polynomial solutions to (0), corresponding to different choices of r(A; fix, . . . , 

Hence the Bethe Ansatz gives C]^_ 1 eigenstates of t(X). On the other hand the Hilbert 
space of the system can be identified with Hi ® • • • ® TCn, where Hi, i = 1, . . . , N 
is a Hilbert space of representation of osp(l,2) corresponding to the z-th particle. All 
eigenstates of i(A) are constructed as linear combinations of basic states in H\®- • -®Hn, 
the latter obtained by the repeated actions of V7's, i — 1, . . . , N on the vacua |0) r The 
number of all possible eigenstates of t(X) at the n-th level equals the number of all possible 
distributions of n numbers between N intervals and hence reads 

This means that still — C^ f _ 1 = C^ -1 states need to be constructed. Using an 
argument similar to the one that lead to [t(X), S z ]- = 0, we can show that [t(X), V + ]- = 0, 
where V+ = E?=i K- The 

remaining states at the n-th level are obtained by the action 
of V+. on each state of the n — 1-th level. 
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Since V + commutes with t(A), it is clear that at the n-th level the full spectrum of 
t(X) consists of 

r (A), r(A;M 1} ), r(A;Ai} B) ,...,ri B) ), 
corresponding to the states 

(v + r\o), (y+^v+oJpm, v + (^\...,^)\o). 

To illustrate, how the theory of orthogonal polynomials can be used to solve the 
osp(l, 2) Gaudin model, we look at the case N = 2, e\ = — e 2 = e, e > 0, s\ = S2 = s. 
Without the loss of generality we can set e = 1. We construct the spectrum of t(X) 
explicitly for this simple model. The Lame equation fl2~j|) can be now written 

(1 - A 2 K(A) - 8s\q' n (X) - 4(r (A) - r(A; p u . . . , /z n ))(l - A 2 )g n (A) = 0, (28) 

where q n (X) is a polynomial of degree n. Eq.(|28|) has a polynomial solution of degree n, 
provided that 

n(n + 8s — 1) 



r(A; m, . . . , fi n ) = r (A) + 



4(a 2 - 1; 



... n(n + 8s - 1) ... , . 

= r (A)+ ^a(A). (29) 

If (P^) is satisfied, then 

g n (A) = p( 4 - 1 ' 4s - 1 )(A), (30) 

where P^ 4s_1 ' 4s_1 ' ) (A) denotes the Jacobi polynomial. The numbers //j, i = 1, . . . , n are 
then zeros of the Jacobi polynomial (|3H|). The full spectrum of t(X) at the n-th excitation 
reads 

{to A + 1 - -a(A); fc = 0, 1, . . . , n}. 

8sA 

The osp(l,2) algebra can be represented in a Fock space of the harmonic oscillator 
as follows |H| 

Sl = \[a t ,ah} + , SL = -\al S\ = \{aU) 2 

n = c ~7i> vl = ~w l = l ^ (31) 
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where aj, a\ are destruction and creation operators, [a*, a' 3 -]_ = 5^ and c 2 = 1, C1C2 = 
— c 2 Ci. In this case s = \ and the Jacobi polynomials ( p0|) become the Legendre polyno- 
mials P„(A), and the spectrum of t(X) reads 

fc(fe + l)(A 2 -l)-l 
*- 4(A 2 - 1) 2 ' ~~ ' ' 

so that the eigenvalues of hamiltonian say, at the n-th level are k(k + l)/8, k = 

0,1,. ..,71. 



V. THE DICKE MODEL 

Here we present another model which satisfies the osp(l, 2) Gaudin algebra and hence 
can be solved in the way described in Section III. This is the model of a bosonic and a 
fermionic oscillator coupled to the osp(l, 2) Gaudin system. Let [b, = 1 be a bosonic 
oscillator and let [f,P]+ = 1 be a fermionic oscillator. It is easy to check that the 
operators 

A(X) = ~X + A°{\), B(X) = b + B°(X), C(A) = -6 f + C°(A) 

V-(\) = ~\f + K°(A), V + (X) = \f + V°(\) (32) 

satisfy the algebra (0) provided that A°(X), B°(X), C°(A), K°(A), V°(\) satisfy it and 
b, f oscillators (anti-)commute with them. We assume that A°(X), B°(X), C°(A), V®(X), 
V+{X) correspond to the iV-particle osp(l, 2) Gaudin model and we show that the Bethe 
Ansatz in this case gives all the eigenstates of t(X). In other words we show that the 
Bethe equations (PT|), which read explicitly in this case 

1 N s k 1 n 1 

-o^ + £ — — = -nE — Z — ' i = l,--.,n, (33) 

k=l k 3=1 r'i 

have precisely C^ +1 = solutions. Using the classical method of Stieltjes (cf. [p^] ) 

we consider the function 

Z(ji l ,...,/j, n ) = e-^ti^ J] l^-ejfcl^ni^r-^l, (34) 

h=l,...,N r<s 

j = l,...,n 

which can be defined in C^ +1 regions in R n corresponding to all possible distributions 
of n variables fa between N +1 intervals, (— 00, e n ], [e n ,e n _i], . . ., [ei,+oo). Let us fix 
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one of these regions. It is clear that at the boundary of this region Z vanishes, and since 
Z > 0, it has a maximum inside the region. At this maximum \i r ^ fi s , if r ^ s, and 
|p: = 0, i = 1, . . . , n, i.e. Eqs. fl33|) hold. Since there are C^ +1 regions, there are Cjy +1 
different maxima, hence (^) has at least C% +1 real solutions. Knowing that there can 
be at most C^ +1 solutions to ([33|) we have also proved that there are precisely C^ +1 
maxima of (0). 

Now we look at the simplest non-trivial case of the Dicke model. We assume that 
iV = 1, ei = 0, and we write osp(l,2) generators in the oscillator representation (j3l[). 
In the purely bosonic case this model corresponds to the system of two oscillators with 
the frequency ratio 2 : 1. In the osp(l, 2) case we have the model of two such oscillators 
with, in addition, an ordinary and one Majorana fermionic oscillators. The generating 
function £(A) takes the form 

t(\) = -\ 2 --H--GX- 1 -—X- 2 , (35) 
v ; 4 2 2 16 v ; 



where 



H = ^[o,ot] + + [ft ) 6t] + _I[/,/t] 



G = ba^ 2 + &V - -^fca^ + ^=fca. (36) 

v 2 v 2 

The use of ( p2|) gives immediately the eigenvalues of t(X) in the form 

r(X; fH, ■ ■ ■ , /in) = ^A 2 - ^A- 2 - \{n + 1) - {\ (37) 

Therefore the hamiltonians H and G have the eigenvalues 

1 n 1 

E n = n + 1, g n =~J2- (38) 

Z i f^i 

respectively. The eigenvalues g n need more discussion. From the symmetries of the 
function Z(ji\, . . . , fi n ) (0), we immediately deduce that if /ij is a solution to (|33|) so is 
—fit, hence the spectrum of G is closed under the change of sign. In addition it contains 
for n even. So the spectrum of G is fully characterised by the positive eigenvalues. On 
the other hand g n can be derived directly from the differential equation (|24]), which now 
takes the form 

Ag^A) - (2A 2 - l)q' n {\) + 2{n\ + g n )q n {\) = 0. (39) 
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Here we have used the explicit form fl38|) of the eigenvalue E n of H . Making the ansatz, 
QnW — J2k=o a k^ k j a n = 1 we can derive the recurrence relations for a^, 

On— 1 Qn 

(n - k) 2 a n _ k + 2{k + 2)a n _ fc _ 2 + 2# n a n _ fc _i = 0, k = 0, 1, . . . , n - 1. (40) 

Solving ( pE0"D for n = 1, 2, we find that the eigenvalues of G are and 0,±V3 respec- 
tively. 

We can represent bosonic oscillators a, b in terms of the canonical coordinates, Xi, Pi, 
i = 1,2, [xi,pj]^ = iSij, as follows 



a = -j^{x 2 + ip 2 ) 



h = xi + -pi. 



We can also represent fermions /, c in terms of the Pauli matrices Oi 

c = cr 3 , /= -fa -ia 2 ). 
In this representation, the hamiltonians ( ftSD take the form 

, f pI + pI + axI + xI + i o \ 



(41) 



(42) 



H 



G 



\ + p\ + 4x? + x\ - 1 / 

/ -2x\p\ + pi [x 2 , p 2 \ + + 2x\x\ -x 2 - ip 2 

\ 



(43) 



-x 2 + ip 2 -2xipl+p 1 [x 2 ,p 2 ] + + 2x 1 xlJ 

It is known that, in the purely bosonic case, Schrodinger equations corresponding 
to G and H can be separated in parabolic coordinates yi, y 2 , 



2xi = y\- y%, 



x 2 = yiy 2 . 



(44) 



It is also known that these are the only coordinates other than Cartesian's, in which the 
bosonic system separates. We can try to implement this procedure in the osp(l, 2) case. 



First we write hamiltonians (43) in parabolic coordinates 
1 ( 

H 



o 



V 



o 



a 2 



(45) 



G 



vlvl 



1 d 2 



dyi yi dyi 



2(y? + yl) [ JL(JL o - y \ + JL a _ 2) 

\ yiy2 K yi oyi 31 j/2 oj/2 W2:/ 
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1 9 2 



4 _ 1 9 2 , „.4 
^ 2 



We look for vector wave functions 



*n (2/1,2/2) = 

which solve the Schrodinger equations 



^(2/1,2/2) 

*n (2/1,2/2) 



) 



TO 
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/I 



9n^ 



(46) 



It can be easily seen that the first of Eqs.(^) can be solved, assuming that ^(2/1,2/2) 
separates, i.e. 



where c± are constants. The equations for $~ are the same as for hence we 

need to consider only equations for say. Eqs.(^8|) are precisely identical with the 
equations one obtains in the purely bosonic case. From this we immediately know, that all 



( f47l) into the second of Eqs.(|46|) we immediately see that if $±, $± satisfy ©, they 
cannot satisfy the Schrodinger equation for G. One has to combine different solutions 
of ((H) to obtain the function *Sf n (yi, 2/2)- Knowing that the parabolic coordinates (|44l) 
are the only ones, for which separation of variables in the first of Eqs. (|46|) is possible, we 
immediately deduce that we cannot separate (^) changing only the position variables. 

One can nevertheless expect that the separation of (f46|) is possible, but a more general 
canonical transformation of the Cartesian variables has to be considered. An indication 
comes from the analysis of the bosonic case. In that case, the zeros of corresponding 
functions $ n can be identified with the zeros of the polynomial q n , satisfying the bosonic 
version of Eq.(p9|). It is therefore reasonable to hope that there are coordinates, in which 
Schrodinger equations (pElf ) give such a natural meaning to Eq.(^). 
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(47) 



(48) 



independent solutions to the first of Eqs. (|4*B]) can have the separated form fl4"7]). Inserting 
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